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Abstract 

We calculate the branching ratios of the — J- n^l'^i'i'y {I = e, /i) decays, and 
the T-odd triple momenta correlations ^ = q- \pi xpVJ/M^, due to the electromagnetic 
final state interaction, in these processes. The contributions on the order of uj^^ and 
LO^ to the corresponding amplitudes are treated exactly. For the branching ratios, 
the corrections on the order of uj are estimated and demonstrated to be small. We 
, compare the results with those of other authors. In some cases our results differ 

Q I considerably from the previous ones. 

<D ' 
Q 

1. The — i- -K^l^vi'-y {I = e, fi) decays were earlier studied theoretically in 
Refs. [H El S]. Therein the branching ratios of these decays were calculated. Besides, 

^ . in Ref. |2] the T-odd triple momenta correlations ^ = q-lpi x 'Pn\/M^ were considered, as 

^! induced by the electromagnetic final state interactions; here and below is the kaon 

Qh| mass, g, p;, are the momenta of 7, , 7r°, respectively. In principle, these triple cor- 

' relations can be used to probe new CP-odd effects beyond the Standard Model, which 

could also contribute to them. 
CS) ■ Here we calculate anew these effects. Our results confirm essentially some previous 

! results and disagree considerably with other ones. 

Ti^ij- \ In the theoretical analysis of radiative effects in the discussed processes, the treatment 

of the accompanying radiation, which gives the effects on the order of uj~^ and uo^ (the last 
ones originate from the radiation due to the lepton magnetic moment), is straightforward 
(here and below u is the photon energy). As to the structure radiation contribution on 
the order of w^, it is also under control, due in fact to the gauge invariance [4J. The 
contributions on the order of cu (and higher) depend directly on the photon field strength 
F^jj (and its derivatives), and cannot be fixed in a model- independent way. We assume 
^ ■ that the corrections on the order of uj and higher are relatively small. And indeed, more 

■ quantitative arguments presented below demonstrate that such contributions into the 

discussed branching ratios do not exceed few percent. 

2. At the tree level, the — ?■ Ti^l^ui'-f decay is described by the Feynman graphs in 

Fig.m 

The matrix elements for diagrams la and lb look as follows: 

Mia = sin6'ceM^7a(l + 75) ^' — —e*vi[f+{t) ■ {pk + Pt,)^ + f-{t) ■ {pk -P7,)a\ 
V 2 '^Pll 

= — =sin6'ceM^7«(l +75) I — ^- Mf+it) " {.PK+P-n)a + f-if) ■ {.PK-Pi.)a\, (1) 
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Figure 1: The tree diagrams 

G PK(i* 

Mib = 7= sin 6'ceu^7„(l + 75)^^/ [/+ (f) ■ {pk - q + Pn)a + f~ (f) ■ (pk -q-P7r)a] ; (2) 

V2 PkQ 

here G is the Fermi couphng constant, 6c is the Cabibbo angle, e is the elementary charge 

(e > 0), t = {pk — P-kY) t' = {pk ~ Q ~ P-kY') the lower indices attached to the matrix 

elements match the corresponding Feynman diagrams. 

Usually the dependence of the form factors /+ and /„ on the momentum transfer t is 

described by formula 

/±W = /±(0)(l + A±^). (3) 

The experimental data are adequately described by Eq. ([2]) with A+ ~ 0.03 for / = /i and 
/ = e; A_ = for / = /i [S]; A_ for / = e is unknown, but one may assume that it is also 
close to zero. 

In the decays the ratio X±t/m'l is small, X±t/m'l < 0.1, so one can put f±{t) = 
/±(0). Since the ratio ^(0) = /_(0)//+(0) ~ 0.1 is also small [6|, one can neglect /_(0) 
with the same accuracy. 

Thus, our expressions ([T]), ([2]) simplify to 

Mia = ^sin6'ce/+(0) ■ (px + P7T)aU^lail + I5) ( — \- 7; — ) (4) 

V2 \Piq 2piqJ 

G PK(i* 

Mu = y= sin6'ce/+(0) ■ {px - q + Pn)a Uuja{l + lb)vi. (5) 

v2 PKq 

However, the sum of diagrams la and lb is not gauge invariant: it does not vanish 
under the substitution e* — )■ q. To restore the gauge invariance, one should add the third 
diagram where a photon is directly emitted from the vertex (see Fig. This contact 
amplitude has no single-particle intermediate states, and therefore is on the order of cu" 
and higher. The contribution ~ 0;°, as derived with the Low technique [1], is 

Q 

Mic = --j= sin 9cef+{0) ■ e*Mi,7„(l + 75)t;/. (6) 

Thus, the model-independent gauge invariant tree amplitude of Kjl^^ decay, including 
only terms on the order of and (but all of them!), is 

Ml = Mia + Mu + Ml, = sin Ocef+iO) ( {pk + p0aM,7„(l + -f5)vi ' ^'"^ ^^"^ 



V2 I \Piq PKq 

+ iPK + P7v)aU,.lail + Ib)^, ^/ + ( -^^qa " 6* ) M,.7a(l + 'y5)vi \ . (7) 

2piq \PKq J J 
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This expression agrees with the corresponding formulas in Ref. [Ij and formula (6) in 
Ref. [2] (if our /+(0) is set to its SU{3) value /+(0) = 1/^2 ). 

It is convenient to present the amplitude ([7]) as a sum of gauge invariant contribu- 
tions. They are the "infrared" term corresponding to the sum of the amplitudes of 
accompanying radiation by the kaon and lepton (independent of the lepton magnetic mo- 
ment), the magnetic term Mmag which is the amplitude of spin-dependent accompanying 
radiation of the lepton magnetic moment, and the Low term M^ow'- 



M, 



IR 



G_ 



PiQ 



M, 



mag 



G 

-j= sin 9cef+{0){pK + Pn)aU,y7a{^ 



Low 



G . . , PRC* 

— =smece/+(0) Qa - e, 

V2 \PKq 



75 k — '^h 
^Piq 



(8) 
(9) 
(10) 



The results of calculation for branching ratios are presented in Table [H here the 
following cuts in the kaon rest frame are used: w ^ 30 MeV and ^ 20°. 





/ = n 


/ = e 


Bijnens et al. \[\ 


1.9 X 10"^ 


2.8 X 10"^ 


Braguta et al. [2J 


2.15 X 10-5 


3.18 X 10"^ 


present work 


1.81 X 10-5 


2.72 X 10-^ 


experimental values 


(2.4 ± 0.5 ± 0.6) X 10-5 [7J 
(1.58 ±0.46 ±0.08) X lO'^ j2| 


(3.06 ±0.09 ±0.14) X 10-^ [9J 



Table 1: Branching ratio of — )■ ir^l^ui'-f decays 



The accuracy of our results can be estimated as follows. The leading corrections to 
them are due to the structure radiation from the hadronic vertex. They are proportional 
to the photon field strength, i.e. are on the order of u. There are good reasons to believe 
that these corrections are less than the Low structure amplitudes which are on the order of 
u^. The Low contributions (including of course their interference with the accompanying 
radiation) to the discussed — 7r°/i"'"z/^7 and — )■ vr^e^z/eT branching ratios, accord- 
ing to our calculations, constitute —0.24 x 10"^ and —0.12 x 10"^, respectively. Thus, 
we estimate the accuracy of our results as ±0.2 x lO^^ and ±0.1 x 10"^, correspondingly. 
Let us note also that corrections to the quoted results derived in Ref. [H |3] in the chiral 
perturbation theory are of similar magnitude. 

As mentioned, additional corrections on the level of 10 - 15% to the branching ratios 
originate from our neglect of the form factor /_(t) and of the t-dependence of f+{t) 
(the last correction is certainly positive). As to the relative accuracy of our numerical 
integration over phase space of final particles, it is about 1%. 

To compare properly our results with those of Refs. jTl [2] one should keep in mind that 
now the experimental values of some quantities are known with better accuracy. Indeed, 
we use sin6'c/+(0) = 0.217/a/2 in our calculation, and, as far as we can see, in Refs. [H |2] 
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the corresponding value is 0.22/a/2. Substitution of one of these values for another alters 
the results by about 3%. 

Thus, our results for the branching ratios agree reasonably well with those of Ref. pLj. 
There is however some disagreement between our results and those of Ref. |2]. 

3. The T-odd triple momenta correlation = q ■ [pi ^ Pt^]/M^ in the — t- n^l^uij 
decays arise from the interference term 2Re(M*742) in the decay rate; here Mi is the tree 
amplitude and A2 is the anti-Hermitian part of the one-loop diagrams presented in Fig. [21 

One can easily demonstrate that A2 is generated only by attaching the intermediate 
photon to diagram la. The on-mass-shell intermediate particles in Fig. |2]are marked by 
crosses. 




b d f 

Figure 2: The one- loop diagrams 



The anti-Hermitian part of the sum of one-loop diagrams is written as 

n 

The Compton amplitude entering this expression is (see Fig. [3]) 

Mfn = + Msb = e VpCk— e vi + e VpC — - — e^fz. (12) 

zpiq —ZpiK 
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-P -P -Pi -P -Pi 

a b 



Figure 3: The Compton scattering diagrams 



As to M*„, it is the same tree amphtude ([7j), up to the change of some notations: 
M4 = Mn^ = Ml I ,^k, = sin e,ef+iO) | (pk + P.)aW.7a(l + ^^^^ 



V2 t VPi'? 



+ + P7r)QMi/7a(l + 75)7r^^P + T^a " ) Mi.7a(l + 75)t;p } ■ (13) 



The sum over n in formula (|TT|) includes the summation over the polarizations of inter- 
mediate particles, and the integral over the phase space with 

dp=^^6^'\p+k-P). (14) 

The details of the calculation of the T-odd correlation = q - [pi ^ p.,r]/Mf^ are given 
in Appendix. Here we wish to note only that our result for ^-odd term |M|^^^ in the 
interference 2Re(M*y42) agrees with formula (12) in Ref. j2] up to the factor 1/a/2 which 
is obviously omitted therein (it follows directly from formulas (6) and (10) of Ref. [2]). 

In fact, what is really measured experimentally, is not the T-odd triple momentum 
correlation C, by itself, but the asymmetry 



/(|Mi| 




|2 + |M|2,,)d$e<o 


ami 







(15) 



/|Mi|2d<l>g>0 



induced by this correlation; here and are the numbers of events with ^ > and 
< 0, and integration is performed over the phase space of the final particles. 

The results for the asymmetry are presented in Table [2J The relative accuracy of 
our numerical integration is about 1%, as it was the case with the branching ratios. 

Here however, as distinct from the problem of branching ratios, the contribution of the 
Low term is large, quite comparable numerically to the contributions of the accompanying 
radiation which are on the order of uj~^ and u^. So, here it is difficult to estimate reliably 
the relative magnitude of the structure radiation contribution proportional to u, i.e. to 
estimate reliably the true accuracy of thus derived results for the asymmetry A^. We note 
here that corrections to the value of the discussed correlation, derived in Ref. [TU] within 
the chiral perturbation theory, are very small. 

Our results for A^ exceed those of Ref. [2j, obtained in the same approximation, by 
a factor of 1.5 - 2. If one included the factor 1/a/2, probably lost in the calculations of 
Ref. |2j (see the remark above), it would make the disagreement even worse. 
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/ = /i 


/ = e 


Braguta et al. [2j 


1.14 X 10"^ 
w > 30 MeV, > 20° 


-0.59 X 10"^ 
w > 30 MeV, > 20° 


present work 


2.38 X 10"^ 
a; ^ 30 MeV, ^ 20° 


-0.93 X 10"^ 
w ^ 30 MeV, 9i.f ^ 20° 


experimental values 


-0.03 ±0.13 [llj 
5 < a; < 30 MeV 


-0.015 ±0.021 [9J 
a; > 10 MeV, 0.6 < cos^,^ < 0.9 


present work 


0.50 X 10-^ 
5 < a; < 30 MeV 


-0.30 X 10"^ 
w > 10 MeV, 0.6 < cos^i^ < 0.9 



Table 2: in -> 7r°/+z//7 decay 



The asymmetry was measured experimentally [HI but cuts imposed therein 
differ from those used in Ref. P]- Therefore, we have calculated also the asymmetry for 
the corresponding kinematical regions (see the last line in Table [2]). 

As distinct from the situation with the branching ratios, all the theoretical results for 
the triple correlations are unfortunately far away from the real experimental sensitivity. 
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Appendix A 

In this section we give the list of the integrals that contribute to the A2: 

dp = ao = ^(l-'-^]; [AA] 



2 







-pi) ! 






ml 















kf,dp = apP^, where op = - f 1 - -^^ j ; (A.2 

/ -\dp = bo = -^In ( ^ ) ; (A.3) 
J Pik 2piq \mfj 

J ^dp = B,^ = hipi^ + bpP^ , (A.4) 
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here bi and bp are the solutions of the system of equations 

birrif + bp{piP) = uq, 
bi{piP) + bpP^^bo{piqy, 



/k k 
-^dp = ^2^^ = b2gu.u + biipi^pi,, + bppP^Pr, + bip{pi^P^ + P^p;,.), (A.5) 

here 62, hi, bpp, and bip are the solutions of the system of equations 

462 + bum^ + bppP^ + 2bip{piP) = 0, 
&2 + burn'i + bipipiP) = 0, 
bpp{piP) + bipm] = ap, 
&2 + bppP"" + bip{piP) = bp{piq)] 



f 1 ^ ^TT 1 ( pKP+^/{pKP?-mlP^ ^ 

J PKk ^ 2 ^(p^P)2 - mlP^ \ pkP - JiPKPY - m%P^ 



(A.6) 



/k 
^^dp = Ci^ = ckPki^ + cpP^, (A. 7) 

here and cp are the solutions of the system of equations 



CKm\ + cp{pkP) = ao, 
ck{pkP) + cpP^ = Co(p;g); 



/ , ,J ,, ^P = ^0 = , ' Inf ^^^^IV^^^M^ ; 

J {Pik){pKk) 2piq JipKPiY - mlm^ \ PkPi - \/{vkViY - mlm^ / 



(A.8) 



/ 



2piq ^J^PkPiY - ■m\m} \PkPi - \/{pkPiY - mj^rn^ 
k 

■j—rrdp = -Di^ dxPRfx + dipi^ + dpP^, (A.9) 



{pik){pKk) 

here di^-, d;, and dp are the solutions of the system of equations 

dRrr^K + di{pKPi) + dpipxP) = bo, 
dxiPRPi) + dim] + dp{piP) = Co, 
dKiPxP) + ciKPi^) + c?P^' = c?o(Pig); 



/ 



k k 

-dp = P>2^,. = d2gij.u + dKKPKpJ)Kv + diiPlixPlu + dppP^P,, 



{Pik){pKk) 

+ dxiiPKuPiu + PiuPkv) + dKp{pKu,Pu + PuPku) + dip{pi^P^ + P^^Pi^), (A.IO) 
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here ^2, dxK, dih dpp, dxi, dxp, and dip are the solutions of the system of equations 



4(^2 + dKKfriK + diim^i + dppP"^ + 2dKi{pKPi) + 2dKp{pKP) + 2dip{piP) = 0, 

d2 + diimf + dxiiPRPi) + dip{piP) = 0, 
dKK{PKPi) + dKim^i + dKp{piP) = ck, 
dppipiP) + dxpiPKPi) + dipin] = cp, 
d2 + dKKm% + dxiiPKPi) + dxpiPRP) = 0, 
dppipxP) + dKpm\ + dip{pKPi) = bp, 
d2 + dppP'^ + dKp{pKP) + dip{piP) = dp{piq). 

All these expressions agree with the analogous ones in Ref. [2]. 

Appendix B 

The A2 can be expressed as follows: 



^2 = ^E^/«^' 



87r2 



+ ^3 



(B.l) 



where 



ie^ G 



^3a-/ij = ^^M3aM7jj|g^fc, = — — sin6'ce/+(0)(pi^ +p^)c,M^7«(l + 75) 



co(P - mi)pK - CiPk H (ao(-P - "^/) - apP) 

Piq 



P — mi ^ ,^ ^ 
-e*vu B.2 

^Piq 



A 



ie^ G 



Sa—maq ^ o 



^ M3„M„„<;| g^fc, = ^ sin ecef+{0){pK + P^)aUulai^ + 75) 



2ap(p2 + 2miP) 



P — rrti ^ ,^ ^ 



A 



3a— Low 



e— >ej; 



sin6'ce/+(0)u^7a(l + 75) 



X 



8712^2 

ao(P - mz)7« - apP'ja - Cia(P - mi)PK 



l^^r.,; (B.4) 



8 



Asb-iR = MsbMml q^k, = -^—E- sin 9cef+{0){pK + P7c)aUu'ya{^ + lb) 



1 

X - 

2 



(P - m^e* ( — ((p; - mi){boP - Bi) - B^P) - do{pi - mi)pK + D,pK 
{Bie*(pi - mi)P - B2^,,i^e*{pi - mi)-^^ - S2Mi^7/^e*7^P) 



+Die*{pi-mi)pK- D2^ui,j,e*-f^PK vi; (B.5) 



i s-^ ie^ G 

Mh-mag = ^ M2,hMmag\ q-^k, = Sin 6lce/+(0)(pK + p^)afii/7a(l + lb) 



1 



X 



B2^ul,j.Pie*lu + 4miS2^^7^e* - Bipie*P - 4miBi(pie*) - mf^ie*] v,; (B.6) 



^3b-Lo«; = -3-^ V AfsfeMLo^lq^fc, = — sin 6'ce/+(0)M,.7a(l + 75) 

OTT^ ^-^ Pl^p, STT^ -v/S 

X - [(P - mi)e*{{pi - mi){DiaPK - &07a) - D2a,j.l^,PK + Pl7a) 

+Bie*{pi - mi)-fa - -B2^t/7/*e*7,.7a] ^^i- (B.7) 

Appendix C 

The interference term 2Re(M*A2) can be expressed as follows: 

2Re(M*A2) = 2Re((M;^ + M;„^ + M^^) 

We keep here only terms odd in ^. Let the \M\jj^_^^_jj^ be the ^-odd terms in 
2Re(M;^^3„_/R), the |M|fR_3„_^„g be the ^odd terms in 2Re(M;^^3„_„„g), etc. 

Therefore, the ^-odd term |M|^j^^ in the interference term 2Re(Mj*742) is 

l-^lodd — l^lm-Sa-IR + \ ll R-3a-mag + \ ^ I'l R-3a- Low + l^lmSb-IR 

+ l^lm-Sb-mag + \ ^ R-3b- Low + \^\mag-3a-IR + \-^\mag-3a-mag 
~^ \-^\mag—3a—Low ~^ \-^\mag—3b—IR 

+ \M\ 

mag— 3b— mag 

+ \M\ 

mag— 3b— Low 

(C.2) 

+ \^\\ow-3a-IR + \^\'Low-3a-mag + \^\\ow-3a-Lmi] + |-^lLo«)-36-7il 
+ \-^\\cyw-3b-mag + |-^lio«)-36-Lo«)' 
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where 



X {3m K -^l- '^{PkPi - PkP-k + Pxq + PiP^ + Ptt^)) 

X ((2ao - ap)m;^ + CKm\{piq) - cpmf{pKq) + 2{cp - Co){pkPi + PKq){Piq)); 

(C.3) 



X ap{mf - 2piq){2,m\ -ml- 2{pkPi - PkPtt + Pxq + PiPn + Pnq)); 

(C.4) 



X ((cx + 2cp){{piqfm\ + {pKqfmf + {pKPi){PKq)m^i) 

+ 2ck{pkPi + PKq){piq)m\ - cp{pkPi + PKq){piq){2piq + m^) 

- 2cK{{pKPif + {PKqf){piq) - 4:CK{pKPi){PKq){piq) 

+ (oo + ap){piq - 2pKq)mf + 2{ap - ao){piqf); 



X (-3m|- + + 2(p7^pi - + + PjPtt + P^q)) 

X (262 + {2bo - 5bi - 5bp + bu + 4bip + 3bpp)m] (C.6) 

+ 2{bo-bi- 3bp + bip + 2bpp){piq) + (rf^ - 2dKp){piq)m\ 

+ (-ciz - 3dp + 2rf/p + 2dpp)mf{pKq) + 2dKp{pKPi){PKq) 

+ 2{-dQ + di + 2dp - dip - dpp){pKPi){piq))] 



Vx/2 



T^^{PKq){piq) 



X (-3m|; + + 2(pi^p/ - pxPn + Pxq + M^tt + P-Kq)) 
X ((6; + 6p - 26;p - 26pp)m;2 + 2(6p - bpp){piq)); 



(C.7) 
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X {-2dKK{piq)mK + '^{dm - 2dpp){piqfm\ 

+ {-dKK + 4rfKi + 4:dKp){pKq)mfm'^j^ + 2{du - dpp){piq)m1m\ 

+ 2{2dK + dKK - 4:dKi - '2,dKp){pKPi){piq)m\ + 2dKK{PKq){piq)m\ 

+ 2(-rfj^; - UKp){pKqfm^i + 4(-26p + fc^p + hpp){piqf 

+ 4(-di - dp + dii + 2dip + dpp){pKPi){piqf 

+ {-dii - Adip - 3dpp){pKq)mf + 2(-6fep + 4bip + 46pp - 3d2){pKq)'mf 

+ 2{dK - ^di + 2dp - dxi - 2dKP + 2dii + Mip + 2dpp){pKPi){pKq)mf 

+ 2(-3fep + 6zz + 26zp + bpp){piq)mf 

+ 4:{-dK + 2di + dxi + dKP - ^du - 2dip){pKPif{piq) 

+ '^h{piq) + 2(-rf, - rfp + dii + 2dzp + dpp){pKPi)iPiq)mf 

+ 4(46p - 26ip - 26pp + d2){pKPi){,Piq) 

+ 2(-(iii - Mip - Mpp){pKq){piq)m] 

+ 4(-rfi^ + rfxi + 2rfxp + '2dpp){pKPi){pKq){piq) 

+ 4(6o - hi){pKq)m] + 2(6o - ?>hi){piq)mf + 4(6o - &z)(Pig - '2pKPi){piq)); 

(C.8) 

l^lL,-3a-7ii= (^^sin^,e/+(0)j -^^(cp(m2-2p,g) + 2co(m)) 

X (3m|- - ml - 2{pkPi - PxPn + Pxq + PzPtt + Ptt?)); 

l^lL.-3a-n.a. = 0; (C.IO) 



l^lL.-3„-.o.= (^sin^.e/,(0)) 

X {2cK{piq)m\ - 2cp{piqf + (cx + 2cp){pkPi + PKq)mf 

- cp{piq)mf - 2ck{pkPi + PKq){Piq) + 4(ap - ao){piq) - 2(ao + ap)mf); 

(C.ll) 



|M|L.-3.-..= (^sin..e/,(0)) ^ 

X (3m|- - - 2(pft:p/ - p^Ptt + Pxq + PzPtt + ^Trg)) 
X {dKKm\ + (dii + 2dip + rfpp - 0?; - rfp)m;^ + 2{dKi + c?kp - dK){pKPi))] 

(C.12) 

l^lL.-36-.na. = 0; (CIS) 
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|M|L,-3.-.«.-(|-<'.e/40)) ^ 

X {{da + 2dip + dpp)mf + 4(26/ + 2bp - 2bii - Akp - 2bpp + d2)mf 

+ dKR^n^m^ + 2{-dK - 2di - 2dp + dKi + dKp){pKPi)mf 

+ 2{dKi + 2dKP + 2dii + Adip + 2d,pp){pKq)mf 

+ 2{dip + dpp){piq)m] - 862 + ^h{m] - 2piq) 

- 2{dKK + '^dxi + 2dKp){piq)m\ + 8(6/ + 2hp - bip - bpp){piq) 

+ A{dK + 2di + 2dp - dxi - dxp - 2d// - 4d/p - 2dpp){pKPi){piq)); 

(C.14) 



K 



T^'^{pKq){piqY 

X Mpiq)mf + 2{ap - 2ao){pKPi)mf + (c/^ - 2cq + 2cp){piqfm\ (C.15) 
- 2cK{pKPi){piq)m\ + 2(co - cp){pKPi){piq){2pKPi + Sp^-g - P/?) 
+ cpm/^(p;^p/)(p/g)); 

X {-2ap){2{pKPi)m1 + (p/g)mf + 2(p/g)' - 4.{pKPi){piq))] 

- 2cp{pKPi){piq) + CK{{piq)m\ - 2{pKPif - 2{pKPi){pKq)))] 



V2 / T^^{PKq){piq) 

X i-2dKKipiq)m% + 2{-dK - 2dp + 2dKP + 2d/p + 2dpp){piqfm\ 
+ 2(-2d/ - 2dp + dm + (^/fp + dii + 2d/p + dpp){piq)m'{mj( 

+ 4:d2{piq)m^K + K^k - dxi - dKp){pKPi){Piq)m]^ + 2dKK{PKq){Piq)m]^ 
+ (26/ + 26p - 36// - 66/p - 3bpp)mf + 4(6/ - 6// - bip){pKq)m'{ 

+ 4(-6o + 6/ + 36p - 6/p - 2bpp){piqf 

+ 4(rfo - - 2rfp + d/p + dpp){pKPi){piqf 

+ 4(26o - 56/ - 56p + 36// + 66,p + 3bpp){pKPi)mf 

+ 2(-26o + 56/ + 66p - 6// - 66/p - 5bpp){piq)mf 

+ S{-do + 2di + 2dp - dii - 2dip - dpp){pKPif{piq) 

+ 2{-di -dp + dii + 2dip + dpp){pKPi){piq)m^ 

+ 8(60 - 26/ - 36p + 6// + 36;p + 2bpp){pKPi){piq) 

+ 2{dp + da + dip){pKq){piq)mf - 462(mf - 2{pkPi) + {piq)) 

+ A{-dK + 2dp + dm - 2dip - 2dpp){pKPi){pKq){piq))] 

(C.18) 
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n'^ipKq)ipiq) 

X {{ki + 2bip + bpp)mf - ^h^m] + 4(6pp - bp){piqf + 8b2{pKPi) 
+ A{bi + bp- 2bu - Abip - 2bpp){pKPi)m1 

+ A{bp - bip - bpp){pKqWi + 2{-bi - bp + 2bip + 2bpp){piq)m^); 

(C.19) 



\/2 7 7r2(p^g) 

X {dKK{PKPi)m\ - {dxi + dKp)mfmj^ + dKi{piq)m\ 
+ (-26, - 2bp + 6,, + 2bip + 6pp)m2 
+ 2(-rfx + dKi + dKp){pKPif + 262 + 2d2{pKPl) 
+ {-du - 2dip - dpp){pKq)m?i 

+ 2dKp{pKPi){PKq) + 2(60 - 6; - 26p + 6,p + bpp){piq) 
+ 2(-(i, - dp + du + 2d,p + dpp){pKPi){piq)- 
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